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Abstract 

The reductions of the Heun equation to the hypergeometric equation by polynomial 
transformations of its independent variable are enumerated and classified. Heun- 
to-hypergeometric reductions are similar to classical hypergeometric identities, but 
the conditions for the existence of a reduction involve features of the Heun equa- 
tion that the hypergeometric equation does not possess; namely, its cross-ratio and 
accessory parameters. The reductions include quadratic and cubic transformations, 
which may be performed only if the singular points of the Heun equation form a 
harmonic or an equianharmonic quadruple, respectively; and several higher-degree 
transformations. This result corrects and extends a theorem in a previous paper, 
which found only the quadratic transformations. [See K. Kuiken, "Heun's equation 
and the hypergeometric equation", SIAM Journal on Mathematical Analysis 10 (3) 
(1979), 655-657.] 

Key words: Heun equation, hypergeometric equation, hypergeometric identity, 
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1 Introduction 



Consider the class of linear second-order differential equations on the Riemann 
sphere CP 1 which are Fuchsian, i.e., have only regular singular points |15j . 
Any such equation with exactly three singular points can be transformed to 
the hypergeometric equation by appropriate changes of the independent and 
dependent variables. Similarly, any such equation with exactly four singular 
points can be transformed to the Heun equation. (See [TU1 Chapter 15]. |22I26] .) 
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Solutions of the Heun equation are much less well understood than hypergeo- 
metric functions |3J. No general integral representation for them is known, for 
instance. Such solutions have recently been used in fluid dynamics |7|25| and 
drift-diffusion theory [S]. They also arise in lattice combinatorics |13|17j . But 
it is difficult to carry out practical computations involving them. An explicit 
solution to the two-point connection problem for the general Heun equation 
is not known |24j . though the corresponding problem for the hypergeometric 
equation has a classical solution. Most work on solutions of the Heun equation 
has focused on special cases, such as the Lame equation |!U|19j . 

Determining which Heun equation solutions are expressible in terms of more 
familiar functions would obviously be useful: it would facilitate the solution 
of the two-point connection problem, and the computation of Heun equation 
monodromies. A significant result in this direction was obtained by Kuiken |18j . 
It is sometimes possible, by performing a quadratic change of the independent 
variable, to reduce the Heun equation to the hypergeometric equation, and 
thereby express its solutions in terms of hypergeometric functions. Kuiken's 
quadratic transformations are not so well known as they should be. The useful 
monograph edited by Ronveaux (221 does not mention them explicitly, though 
it lists Ref. |^ in its bibliography. One of Kuiken's transformations was re- 
cently rediscovered by Ivanov [T^j, in a disguised form. 

Unfortunately, the main theorem of Ref. ^H] is incorrect. The theorem asserts 
that a reduction to the hypergeometric equation, by a rational change of the 
independent variable, is possible only if the singular points of the Heun equa- 
tion form a harmonic quadruple in the sense of projective geometry; in which 
case the change of variables must be quadratic. In this paper, we show that 
there are many alternatives. A reduction may also be possible if the singular 
points form an equianharmonic quadruple, with the change of variables be- 
ing cubic. Additional singular point configurations permit changes of variable 
of degrees 3, 4, 5, and 6. Our main theorem (Theorem I3.1j) and its corollar- 
ies classify all such reductions, up to affine automorphisms of the Heun and 
hypergeometric equations. It replaces the theorem of Ref. [18J. 

It follows from Theorem 13.11 that in a suitably defined 'nontrivial' case, the 
local Heun function HI can be reduced to the Gauss hypergeometric func- 
tion 2 Fi by a formula of the type Hl(t) = 2 Fi(R(t)) only if the pair (d, q/a/3), 
computed from the parameters of HI, takes one of exactly 23 values. These 
are listed in Theorem 13.71 A representative list of reductions is given in Theo- 
rem |3~H1 These theorems should be of interest to special function theorists and 
applied mathematicians. We were led to our correction and expansion of the 
theorem of Ref. [THj by a discovery of Clarkson and Olver jB] : an unexpected 
reduction of the Weierstrass form of the equianharmonic Lame equation to 
the hypergeometric equation. In §0 we explain how this is a special case of 
the cubic Heun-to-hypergeometric reduction. 
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The new reductions are similar to classical hypergeometric transformations. 
(See Chapter 3]; also |TDJ Chapter 2].) But reducing the Heun equation to 
the hypergeometric equation is more difficult than transforming the hypergeo- 
metric equation to itself, since conditions involving its singular point location 
parameter and accessory parameter, as well as its exponent parameters, must 
be satisfied. Actually, the reductions classified in this paper are of a somewhat 
restricted type, since unlike many classical hypergeometric transformations, 
they involve no change of the dependent variable. A classification of reductions 
of the more general type is possible, but is best phrased in algebraic-geometric 
terms, as a classification of certain branched covers of the Riemann sphere by 
itself. A further extension would allow the transformation of the independent 
variable to be algebraic rather than polynomial or rational, since at least one 
algebraic Heun-to- hypergeometric reduction is known to exist Extended 
classification schemes are deferred to one or more further papers. 



2 Preliminaries 



2.1 The Equations 



The Gauss hypergeometric equation is 

a + b — c + 1 





ft) 



where a, b, c <G C are parameters. It and its solution space are specified by the 
Riemann P-symbol 



P < 




(2.1) 



where each column, except the last, refers to a regular singular point. The first 
entry is its location, and the final two are the characteristic exponents of the 
solutions there. The exponents at each singular point are obtained by solving 
an indicial equation [T^]. In general, each finite singular point z has ( as an 
exponent if and only if the equation has a local (Frobenius) solution of the 
form (z—Zoph(z) in a neighborhood of z — Zq, where h is analytic and nonzero 
at z = zq. If the exponents at z = zq differ by an integer, this statement must 
be modified: the solution corresponding to the smaller exponent may have a 
logarithmic singularity at z = zq. The definition extends in a straightforward 
way to Zq = oo, and also to ordinary points, each of which has exponents 0, 1. 
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There are 2x3 = 6 local solutions of (fj) in all: two per singular point. If 
c is not a nonpositive integer, the solution at z = belonging to the exponent 
zero will be analytic. When normalized to unity at z = 0, it will be the Gauss 
hypergeometric function 2-^1(0, b; c; z) [TO]. This is the sum of a hypergeometric 
series, which converges in a neighborhood of z — 0. In general, 2-^1(0, c; 2) 
is not defined when c is a nonpositive integer. 



The Heun equation is usually written in the form 

■' " 5 e \ d« a(3t — q 




t(t-l)(t-d) 



u = 0. (S3) 



Here d G C, the location of the fourth singular point, is a parameter (d 7^ 0, 1), 
and a, /3, 7, 5, e G C are exponent-related parameters. The P-symbol is 



P < 




(2.2) 



This does not uniquely specify the equation and its solutions, since it omits 
the accessory parameter q G C. The exponents are constrained by 

a + /3-7-5-e + l = 0. (2.3) 

This is a special case of Fuchs's relation, according to which the sum of the 
2n characteristic exponents of any second-order Fuchsian equation on CP 1 
with n singular points must equal n — 2 |2"T] . 

There are 2x4 = 8 local solutions of (ft) in all: two per singular point. If 
7 is not a nonpositive integer, the solution at t = belonging to the exponent 
zero will be analytic. When normalized to unity at t — 0, it is called the 
local Heun function, and is denoted Hl(d,q;a, /3,j,5;t) j22|- It is the sum of 
a Heun series, which converges in a neighborhood of t = |22|26j . In general, 
Hl(d, q; a, (3, 7, 5; t) is not defined when 7 is a nonpositive integer. 

If e = and q = a/3d, the Heun equation loses a singular point and becomes a 
hypergeometric equation. Similar losses occur if S — 0, q — a/3, or 7 = 0, q = 0. 
This paper will exclude the case when the Heun equation has fewer than four 
singular points, since reducing (fj) to itself is a separate problem, leading to 
the classical hypergeometric transformations. The following case, in which the 
solution of (H3) can be reduced to quadratures, will be initially excluded. 

Definition 2.1. If a(3 = and q = 0, the Heun equation (H^) is said to be 
trivial. Triviality implies that one of the exponents at t = 00 is zero (i.e., 
a (3 = 0), and is implied by absence of the singular point at t = 00 (i.e., 

a(3 = 0, a + (3 = 1, q = 0). 
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The transformation to (Sj) or (f)) of a linear second-order Fuchsian differential 
equation with singular points at t = 0, 1, d, oo (resp. z = 0,1, oo), and with ar- 
bitrary exponents, is accomplished by certain linear changes of the dependent 
variable, called F-homotopies. (See [TU] and (221 §A2 and Addendum, §1.8].) 
If an equation with singular points at t — 0, 1, d, oo has dependent variable u, 
carrying out the substitution u(t) = t~ p (t — l)~ a (t — d)~ T u(t) will convert the 
equation to a new one, with the exponents at t — 0, l,d reduced by p,u,r 
respectively, and those at t = oo increased by p + a + r. By this technique, 
one exponent at each finite singular point can be shifted to zero. 

In fact, the Heun equation has a group of F-homotopic automorphisms iso- 
morphic to (Z 2 ) 3 , since at each of t — 0, 1, d, the exponents 0, ( can be shifted 
to — C,0, i.e., to 0, — C- Similarly, the hypergeometric equation has a group 
of F-homotopic automorphisms isomorphic to (Z2) 2 . These groups act on the 
6 and 3-dimensional parameter spaces, respectively. For example, one of the 
latter actions is (a, b; c) 1— > (c — a, c — b; c), which is induced by an F-homotopy 
at z — 1. From this F-homotopy follows Euler's transformation [21 §2.2] 

2 F 1 (a, b; c; z) = (1 - z) c - a -\F 1 (c ~a,c-b;c; z), (2.4) 

which holds because 2-P1 is a local solution at z — 0, rather than at z — 1. 

If the singular points of the differential equation are arbitrarily placed, trans- 
forming it to the Heun or hypergeometric equation will require a Mobius (i.e., 
projective linear or homographic) transformation, which repositions the sin- 
gular points to the standard locations. A unique Mobius transformation maps 
any three distinct points in CP 1 to any other three; but the same is not true 
of four points, which is why (fj) has the singular point d as a free parameter. 



2.2 The Cross-Ratio 



The characterization of Heun equations that can be reduced to the hyper- 
geometric equation will employ the cross-ratio orbit of {0, 1, d, 00}, defined as 
follows. If A, B,C,D E CP 1 are distinct, their cross-ratio is 



{A B; ^' D) ~ f D -^{c -B) E U °' °° } ' (2 ' 5) 

which is invariant under Mobius transformations. Permuting A,B,C,D yields 
an action of the symmetric group S4 on CP 1 \ {0, l,oo}. The cross-ratio is 
invariant under interchange of A, B and C, D, and also under simultaneous 
interchange of the two points in each pair. So each orbit contains no more 
than 4!/4 = 6 cross-ratios. The possible actions of 5*4 on s £ CP 1 \ {0, 1, 00} 
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are generated by s i— > 1 — s and s i— > 1/s, and the orbit of s comprises 



s, 1/(1 -s), s/(s-l), (s-l)/s, (2.6) 

which may not be distinct. This is called the cross-ratio orbit of s; or, if s = 
(A, B;C, D), the cross-ratio orbit of the unordered set {A, B,C, D} C CP 1 . 
Two sets of distinct points {Ai, Bi, Ci, Di} (i = 1, 2) have the same cross-ratio 
orbit iff they are related by a Mobius transformation. 

Cross-ratio orbits generically contain six values, but there are two exceptional 
orbits: one with three and one with two. If (A, B; C, D) = — 1, the cross-ratio 
orbit of {A, B, C, D} will be {-1, \, 2}. The value -1 for (A, B; C, D) defines a 
so-called harmonic configuration: A, B and C, D are said to be harmonic pairs. 
More generally, if {A, B, C, D} has cross-ratio orbit { — 1, |, 2}, it is said to be 
a harmonic quadruple. It is easy to see that if C = oo and A, B, D are distinct 
finite points, then A, B and C, D will be harmonic pairs iff D is the midpoint 
of the line segment AB. In consequence, {A, B, oo, D} will be a harmonic 
quadruple iff {A, B, D} C C comprises collinear, equally spaced points. So, 
{A, B, C, D} C CP 1 will be a harmonic quadruple iff it can be mapped by a 
Mobius transformation to a set consisting of three equally spaced finite points 
and the point at infinity; equivalently, to the vertices of a square in C. 

The cross-ratio orbit containing exactly two values is {| ± i-^}- Any set 
{A, B, C, D} with this as cross-ratio orbit is said to be an equianharmonic 
quadruple. {A, B, oo, D} will be an equianharmonic quadruple iff A, B, D are 
the vertices of an equilateral triangle in C. If CP 1 is interpreted as a sphere via 
the usual stereographic projection, then by an afline transformation (a special 
Mobius transformation), this situation reduces to the case when A, B,oo,D 
are the vertices of a regular tetrahedron. So, {A, B, C, D} C CP 1 will be an 
equianharmonic quadruple iff it can be mapped by a Mobius transformation 
to the vertices of a regular tetrahedron in CP 1 . 

Cross-ratio orbits are of two sorts: real orbits such as the harmonic orbit, 
and non-real orbits such as the equianharmonic orbit. All values in a real 
orbit are real, and in a non-real orbit, all have a nonzero imaginary part. So, 
{A, B, C, D} will have a specified real orbit as its cross-ratio orbit iff it can 
be mapped by a Mobius transformation to a set consisting of three specified 
collinear points in C and the point at infinity; equivalently, to the vertices 
of a specified quadrangle (generically irregular) in C. Similarly, it will have a 
specified non-real orbit as its cross-ratio orbit iff it can be mapped to a set 
consisting of a specified triangle in C and the point at infinity; equivalently, 
to the vertices of a specified tetrahedron (generically irregular) in CP 1 . 

The cross-ratio orbit of {0, l,d, oo} will be the harmonic orbit iff d equals 
— 1, |, or 2, and the equianharmonic orbit iff d equals | ± i^. In contrast, 
it will be a specified generic orbit iff d takes one of six orbit-specific values. 
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The cross-ratio orbit of {0, l,d, 00} being a specified orbit is equivalent to 
its being the same as the cross-ratio orbit of some specified quadruple of the 
form {0,1, D, 00}, i.e., to there being a Mobius transformation that maps 
{0, 1, d, 00} onto {0, 1, D, 00}. The possibilities are 

d = D, 1-D, 1/D, 1/(1 -£>), D/(D-1), (D-l)/D. (2.7) 

By examination, this is equivalent to {0, 1, d} being mapped onto {0, 1, D} 
by some affine transformation, i.e., to AOlrf being similar to A01.D. The 
corresponding affine transformations t 1— > Ai(t) are 

A 1 (t)=t, l-t, Dt, (D-l)t+l, (1-D)t + D, D(l-t). (2.8) 

This interpretation gives a geometric significance to the possible values of d. 

For the equianharmonic orbit, in which the six values degenerate to two, the 
triangle A01-D may be taken to be any equilateral triangle. For any real orbit, 
AO ID must be taken to be degenerate, with collinear vertices. For example, 
the cross-ratio orbit of {0, 1, d, 00} being the harmonic orbit is equivalent to 
A01<i being similar to A012, i.e., to {0, 1, d} consisting of three equally spaced 
collinear points. This will be the case iff d equals —1, |, or 2, in agreement 
with the definition of a harmonic quadruple. 

2.3 Automorphisms 

According to the theory of the Riemann P-function, any Mobius transforma- 
tion M of the independent variable will preserve characteristic exponents. For 
the hypergeometric equation (jQ), this implies that if M is one of the 3! Mobius 
transformations that permute the singular points z = 0, 1, 00, the exponents of 
the transformed equation (f)) at its singular points M(0), M(l), M(oo) will be 
those of (0) at 0, 1, 00. But if M is not affine, i.e., M(oo) 7^ 00, then (()) will 
not in general be a hypergeometric equation, since its exponents at M(oo) 
may both be nonzero. To convert (I)) to a hypergeometric equation, the per- 
mutation must in this case be followed by an F-homotopic transformation of 
the form y(z) = [z — M(oo)]~ a y(z) or y(z) = [z — M (oo)]~ b y(z). 

Definition 2.2. Aut(fj), the automorphism group of the hypergeometric equa- 
tion (j^J, is the group of changes of variable (Mobius of the independent 
variable, linear of the dependent) which leave (jf)]) invariant, up to parame- 
ter changes. Similarly, Aut(fj) is the automorphism group of (fj). 

Aut(()) acts on the 3-dimensional parameter space of (0). It contains the sym- 
metric group S3 of permutations of singular points as a subgroup, and the 
group (Z 2 ) 2 of F-homotopies as a normal subgroup. So Aut(f)) ~ (Z 2 ) 2 x S3, 
a semidirect product. It is isomorphic to S 4 , the octahedral group j^j. 
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Definition 2.3. Within Aut(fj), the Mobius automorphism subgroup is the 
group -M(()) = f {1} x S 3 , which permutes the singular points z — 0, 1, oo. The 
subgroup of affine automorphisms is A(\)) = f {1} x S2, which permutes the 
finite singular points z = 0, 1, and fixes 00. (It is generated by the involution 
z 1— ► 1 — z.) The F-homotopic automorphism subgroup is (Z 2 ) 2 x {1}. 

The action of Aut(fj) on the 2x3 = 6 local solutions is as follows. |Aut(f))| = 
2 2 x 3! = 24, and applying the transformations in Aut(fj) to any single local 
solution yields 24 solutions of (J^J. Applying them to 2F1, for instance, yields 
the 24 series solutions of Kummer ||. However, the 24 solutions split into six 
sets of four, since for each singular point z e {0, 1, 00} there is a subgroup 
of Aut(f)) of order 2 1 x 2! = 4, each element of which fixes z = z and performs 
no F-homotopy there; so it leaves each local solution at z = z invariant. 

For example, the four transformations in the subgroup associated to z = 
yield four equivalent expressions for 2-^1(0, b; c; z); one of which is 2-^1(0, b; c; z) 
itself, and another of which appears above in ()2.4|) . The remaining two are 
(1 - z)~ a 2 F 1 (a, c-b;c; z/(z - 1)) and (1 - z)~ h 2 F 1 (6, c-a;c; z/(z - 1)). The 
five remaining sets of four are expressions for the five remaining local solutions. 
One that will play a role is the 'second' local solution at z = 0, which belongs 
to the exponent 1 — c. One of the four expressions for it, in terms of 2-F1, is JU] 

2Fi(a, b- c; z) = z^^F^a - c + 1, b - c + 1; 2 - c; z). (2.9) 
2Fi(a, b; c; z) is defined if c 7^ 2, 3, 4, . . . The second local solution must be 
specified differently if c = 2, 3, 4, . . . , and also if c = 1, since in that case, 2-F1 
reduces to 2F1. (Cf. |TJ § 15.5].) When 2 F\ is defined, it may be given a unique 
meaning by choosing the principal branch of z l ~ c . 

The automorphism group of the Heun equation is slightly more complicated. 
There are 4! Mobius transformations M that map the singular points t = 
0, 1, d, 00 onto t = 0,1, d', 00, for some d! E CP 1 \{0, 1, 00}. The possible d' con- 
stitute the cross-ratio orbit of {0, 1, d, 00}. Of these 4! transformations, 3! fix 
t = 00, i.e., are affine. All values d' are obtained by affine transformations, i.e., 
a mapping is possible iff AOlci is similar to A01(f . (Cf. the discussion in § 12.21 ) 
If M is not affine, it must be followed by an F-homotopic transformation of 
the form u(t) = [t - M(oo)]- a u(t) or u(t) = [t - M (oo)]" /3 u(t). 

Aut(f)) acts on the 6-dimensional parameter space of (J^j). It contains the 
group S4 of singular point permutations as a subgroup, and the group (Z 2 ) 3 
of F-homotopies as a normal subgroup. So Aut(5}) ~ (Z 2 ) 3 x 54. It turns out 
to be isomorphic to the Coxeter group V4 [T2"] . 

Definition 2.4. Within Aut(fj), the Mobius automorphism subgroup is the 
group M{fj) = {1} x S A , which maps between sets of singular points of the 
form {0, 1, d' , 00}. The subgroup of affine automorphisms is A(S)) == {1} x S s , 
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which maps between sets of finite singular points of the form {0, l,d'}, and 
fixes oo. The F-homotopic automorphism subgroup is (Z 2 ) 3 x {1}. 

The action of Aut(^) on the 2x4 = 8 local solutions is as follows. |Aut(i})| = 
2 3 x 4! = 192, and applying the transformations in Aut(ij) to any single local 
solution yields 192 solutions of (^J). However, the 192 solutions split into eight 
sets of 24, since for each singular point t G {0, 1, d, oo} there is a subgroup 
of Aut(fj) of order 2 2 x 3! = 24, each element of which fixes t = £q and performs 
no F-homotopy there; so it leaves each local solution at t = to invariant. This 
statement must be interpreted with care: selecting t = d selects not a single 
singular point, but rather a cross-ratio orbit. 

The 24 transformations in the subgroup associated to to = yield 23 equivalent 
expressions for Hl(d, q; a, (3, 7, 5; t), one of which, the only one with no F- 
homotopic prefactor, appears on the right in the identity [22 26j 

HI (d, q; a, j3, 7, 5; t) = HI (1/d, q/d; a, f3, 7, a + /3 — 7 — 8+1; t/d) . 

(2.10) 

(The two sides are defined if 7 is not a nonpositive integer.) The remaining 
seven sets of 24 are expressions for the remaining seven local solutions. One 
that will play a role is the 'second' solution at t = 0, which belongs to the 
exponent 1 — 7. One of the 24 expressions for it, in terms of HI, is 

Hl(d, q; a, 0, 7, 5; t) = t 1 " 7 Hl(d, q; a-7+l, /3-7+1, 2-j, 5; t), (2.11) 

where the transformed accessory parameter q equals q + (1 — j)(e + d5). The 
quantity Hl(d, q; a, [3, 7, 5; t) is defined if 7 7^ 2, 3, 4, . . . The second local so- 
lution must be^ specified differently if jj= 2, 3, 4, ... , and also if 7 = 1, since 
in that case, HI reduces to HI. When HI is defined, it may be given a unique 
meaning by choosing the principal branch of t 1-7 . 

In general, transformations in Aut(i^) will alter not merely d and the exponent 
parameters, but also the accessory parameter q. This is illustrated by ()2.10j) 
and ()2.11|) . The general transformation law of q is rather complicated. Partly 
for this reason, no satisfactory list of the 192 solutions has appeared in print. 
The original paper of Heun ^3] tabulates 48 of the 192, but omits the value 
of q in each. His table also unfortunately contains numerous misprints and 
cannot be used in practical applications fI5\ § 6.3]. Incidentally, one sometimes 
encounters a statement that there are only 96 distinct solutions j4IHI22j . This 
is true only if one uses (|2.1(J|) to identify the 192 solutions in pairs. 
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3 Polynomial Heun-to-Hypergeometric Reductions 

We now state and prove Theorem our corrected and expanded version of 
the theorem of Kuiken ^H] • 

The theorem will characterize when a homomorphism of rational substitution 
type from the Heun equation (|^J to the hypergeometric equation (jQ) exists. 
It will list the possible substitutions, up to affine automorphisms of the two 
equations. It is really a characterization of the .A -orbits that can be mapped 
by homomorphisms of this type to ,4(f))-orbits. The possible substitutions, it 
turns out, are all polynomial. 

For ease of understanding, the characterization will be concrete: it will state 
that AOld must be similar to one of five specified triangles of the form A01-D. 
By the remarks in § 12. 2[ similarity occurs iff d belongs to the cross-ratio orbit 
of D, i.e., iff D can be generated from d by repeated application of d i— > 1 — d 
and d i— > 1/d. The two exceptional cross-ratio orbits, namely { — 1, |,2} (har- 
monic) and {|±i^} (equianharmonic) , will play a prominent role. It is worth 
noting that if KeD = |, the orbit of D is closed under complex conjugation. 

For each value of D, the polynomial map from CP 1 3 t to CP 1 3 z, which 
is denoted R, will be given explicitly when d = D. If is any other value on 
the cross-ratio orbit of D, as listed in (J2.7j) . the polynomial map would be 
computed by composing with the corresponding affine transformation A\ of C 
that takes AOlc? to A01-D; which is listed in (|2.8|) . So if d ^ D, statements 
in the theorem dealing with singular points, characteristic exponents, and 
the accessory parameter must be altered. For example, case 1 of the theorem 
refers to a distinguished singular point do, the mandatory value of which is 
given when d = D. If d ^ D, its mandatory value would be computed as 
the preimage of that point under A\. Similarly, a statement "the exponents 
of t — must be 0, 1/2", valid when d = D, would be interpreted if d ^ D 
as "the exponents of the preimage of t — under A\ must be 0, 1/2". And a 
statement that q/afl must take some value would be interpreted if d ^ D as 
a statement that q/af3 must equal the preimage of that value under A x . 

In the statement of the theorem and what follows, S d = 1 — R. 

Theorem 3.1. A Heun equation (^J) ; which has four singular points and is 
nontrivial {i.e., a(3 ^ or g ^ 0), can be transformed to the hypergeometric 
equation ((}) by a rational substitution z = R(t) if and only if R is a polyno- 
mial, a/3 0, and one of the following two conditions is satisfied. 

(1) AOld is similar to A01.D, for one of the values of D listed in subcases 
la-lc; each of which is real, so the triangle must be degenerate. Also, the 
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normalized accessory parameter q/a/3 must equal one of 0,1, d, which 
may be denoted d . Each subcase lists the value of d when d = D. 
(2) AOlrf is similar to A01D, for one of the values of D listed in subcases 2a- 
2d; each of which is non-real and has real part equal to \, so the triangle 
must be isosceles. Each subcase lists the value of q/afl when d = D. 

Besides specifying D and the value of q/af3 when d = D, each subcase imposes 
restrictions on the characteristic exponent parameters at the singular points 
0,1, d. The subcases of the 'real' case 1 are the following. 

(la) [Harmonic (equally spaced collinear points) case] D = 2. Suppose d = D. 
Then do must equal 1, and t — 0, d must have the same characteristic 
exponents, i.e., 7 = e. In general, either R or S will be the degree-2 
polynomial t(2 — t), which maps t = 0,d to z = and t = 1 to z — 1 (with 
double multiplicity). There are special circumstances in which R may be 
quartic, which are listed separately, as subcase lc. 

(lb) D — 4. Suppose d = D. Then do must equal 1, the point t = 1 must 
have characteristic exponents that are double those oft = d, i.e., 1 — 5 — 
2(1 — e), and t = must have exponents 0, 1/2, i.e., 7 = 1/2. Either R 
or S will be the degree-3 polynomial (t — 1) 2 (1 — 1/4), which maps t = 
to z — 1 and t — l,d to z — (the former with double multiplicity) . 

(lc) [Special harmonic case.] D = 2. Suppose d = D. Then do must equal 1, 
and t — 0, d must have the same characteristic exponents, i.e., 7 = e. 
Moreover, the exponents of t = 1 must be twice those of t = 0,d, i.e., 
1 — 5 = 2(1 — 7) — 2(1 — e) . Either R or S will be the degree-4 polynomial 
4[i(2 — t) — \} 2 , which maps t — 0, 1, d to z — 1 (t — 1 with double 
multiplicity) . 

The subcases of the 'non-real' case 2 are the following. 

(2a) [Equianharmonic (equilateral triangle) case.} D — \ + i^. q/a(3 must 
equal the mean of 0,1, d, and t = 0,1, d must have the same charac- 
teristic exponents, i.e., 7 = 5 = e. Suppose d = D. Then q/af3 must 
equal \ + i^. In general, either R or S will be the degree-3 polynomial 

1 — t/(^ + i^) ' , which maps t = 0, 1, d to z = 1 and t = q/a/3 to z = 
(with triple multiplicity) . There are special circumstances in which R may 
be sextic, which are listed separately, as subcase 2d. 
(2b) D = I + Suppose d = D. Then q/a/3 must equal | + t — d must 
have characteristic exponents 0, 1/3, i.e., e = 2/3, andt = 0, 1 must have 
exponents 0,1/2, i.e., 7 = 5= 1/2. Either R or S will be the degree-4 

polynomial 1 — £/(| + i^p) 1 — + i"^) > which maps t = d, q/a/3 
to z = (the latter with triple multiplicity) and t — 0, 1 to z — 1. 
(2c) D — I + i 11 ^^ . Suppose d = D. Then q/a/3 must equal \ + i^-, t = d 
must have characteristic exponents 0,1/2, i.e., e = 1/2, and t — 0, 1 
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must have exponents 0, 1/3, i.e., 7 = 5 = 2/3. Either R or S will be the 
degree-h polynomial At(t — 1) t — (| + i ^jp) , which mapst = 0, l,q/af3 
to z = (£/ie Zast u>rf/i tnp/e multiplicity) . The factor A is chosen so that 
it maps t = d to z = 1, as well; explicitly, A = —i 202 ^^ . 

(2d) [Special equianharmonic case.] D = | q/af3 must equal the mean 

of 0, l,d, and t = 0, l,d must have characteristic exponents 0, 1/3, i.e., 
7 = 5 = e = 2/3. Suppose d = D. Then q/a/3 must equal | + if. Either 

R or S will be the degree-Q polynomial 4 {[1 — t/(~ + if)] 3 — f } , which 
maps t = 0, 1, d, g/a/5 to z = 1 (t/ie /ast inp/e multiplicity) . 

Remark 3.1.1. The origin of the special harmonic and equianharmonic sub- 
cases is easy to understand. In subcase lc, t ^ R(t) or S(t) is the composition 
of the quadratic map of subcase la with the map z i— > 4(z — |) 2 . In subcase 2d, 
t i — > or S(t) is similarly the composition of the cubic map of subcase 2a 
with z i — > 4(z — |) 2 . In both lc and 2d, the further restrictions on exponents 
make possible the additional quadratic transformation of z, which transforms 
the hypergeometric equation into itself (see [21 §3.1] and [TUj). 

Remark 3.1.2. R is determined uniquely by the choices enumerated in the 
theorem. There is a choice of subcase, a choice of d from the cross-ratio orbit 
of D, and a binary choice between R and S. The final two choices amount 
to choosing affine maps A\ G A{ff) and A 2 G A(i)), i.e., A 2 (z) = z or 1 — z, 
which precede and follow a canonical substitution. 

In the harmonic case la, in which the .A(fj)-orbit includes three values of d, 
there are accordingly 3x2 = 6 possibilities for R; namely, 

R = t 2 , l-t 2 - (2t — l) 2 , 1 — (2t - l) 2 ; t(2-t), l-t(2-t), (3.1) 

corresponding to d = —1, —1; |, |; 2, 2, respectively. These are the quadratic 
transformations of Kuiken [IS]. In the equianharmonic case 2a, in which the 
orbit includes only two values of d, there are 2x2 = 4 possibilities; namely, 

R=[l- t/{\ ± if)] 3 , 1 - [1 - t/{\ ± if)] 3 , (3.2) 

corresponding to d = | ± if. The remaining subcases, with the exception of 
lc and 2d, correspond to generic cross-ratio orbits: each value of D specifies 
six values of d. In each of those subcases, there are 6 x 2 = 12 possibilities. So 
in all, there are 56 possibilities for R. 

Remark 3.1.3. The characteristic exponents of the singular points z = 0, 1, oo 
of (jfjj can be computed from those of the singular points t — 0, 1, d, oo of (1^1) , 
together with the formula for R. The computation relies on Proposition 13.31 
below, which may be summarized thus. If t = to is not a critical point of the 
map t i — ► z — R{t), then the exponents of z — R(to) will be the same as 
those of t . If, on the other hand, t = to is mapped to z = z = R(to) with 
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multiplicity k > 1, i.e., t = to is a k — 1-fold critical point of R and z = zq is 
a critical value, then the exponents of Zq will be 1/k times those of to- 

For example, in the harmonic case la, the map t i— > z takes two of t = 

0. 1, d to either z — or z — 1, and by examination, the coalesced point 
is not a critical value of the map; so the characteristic exponents of those 
two points are preserved, and must therefore be the same, as stated in the 
theorem. On the other hand, the characteristic exponents of the third point 
of the three, t = do, are necessarily halved when it is mapped to z — 1 
or z = 0, since by examination, R always has a simple critical point at t = do, 

1. e., z ~ const + C(t — do) 2 for some nonzero C. (These statements follow 
by considering the canonical d = D case.) So if So denotes the parameter 
(out of 7, 5, e) corresponding to t = do, the characteristic exponents of z = 1 
or z = will be 0, (1 — 5 )/2. R, being a quadratic polynomial, also has a 
simple critical point at t = oo, so the characteristic exponents of z = oo are 
one-half those of t = oo, i.e., a/2, (3/2. It follows that in the harmonic case, 
the Gauss parameters (a, b; c) of the resulting hypergeometric equation will be 
(a/2, /3/2; (5 + l)/2) or (a/2, (3/2; ( a + f3-5 + l)/2). 

In the equianharmonic case 2a, the map t h-> z takes t = 0, l,d to either 
2 = or z = 1; and by examination, the coalesced point is not a critical 
value of the map; so the characteristic exponents of those three points are 
preserved, and must therefore be the same, as stated in the theorem. On the 
other hand, at t = q/a(3, which is mapped to z = 1 or z = 0, R has, by 
examination, a double critical point, i.e., z ~ const + C(t — q/a(3) 3 for some 
nonzero C. So the characteristic exponents of z = 1 or z = 0, since t = q/a(3 
is an ordinary point of the Heun equation and effectively has characteristic 
exponents 0,1, are 0,1/3. R, being a cubic polynomial, also has a double 
critical point at t = oo, so the characteristic exponents of z = oo are one- 
third those of t = oo, i.e., a/3, (3/3. It follows that in the equianharmonic 
case, the parameters (a, b; c) of the resulting hypergeometric equation will be 
(a/3, (3/3; 2/3) or (a/3, (3/3; (a + (3+ l)/3). 

Definition 3.2. A rational map R : CP 1 — > CP 1 is said to map the character- 
istic exponents of the Heun equation J^J) to the characteristic exponents of the 
hypergeometric equation (fj) if, for all to G CP 1 , the exponents of t = to ac- 
cording to the Heun equation, divided by the multiplicity of to t— > z == R(to), 
equal the exponents of z = z according to the hypergeometric equation. 

For example, if to and zq are both finite, this says that if z ~ zq + C(t — to) k to 
leading order, for some nonzero C, then the exponents of z = zq must be those 
of t = to, divided by k. If t = t is an ordinary point of the Heun equation, 
then the exponents of z = z will be 0, 1/k. This implies that if k > 1, z = Zq 
must be one of the three singular points of the hypergeometric equation. 
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Proposition 3.3. A Heun equation of the form can be reduced to a hyper- 
geometric equation of the form (j^J by a rational substitution z = R(t) of its 
independent variable only if R maps exponents to exponents. 

Proposition 13.31 which was already used in Remark 13.1.31 above, is a spe- 
cial case of a basic fact in the theory of the Riemann P-function: if a ratio- 
nal change of the independent variable transforms one Fuchsian equation to 
another, then the characteristic exponents are transformed multiplicatively. 
It can be proved by examining the effects of the change of variable on each 
local (Frobenius) solution. It also follows immediately from Lemma E31 below. 

This lemma begins the study of sufficient conditions for the existence of a 
Heun-to-hypergeometric transformation. Finding them requires care, since an 
accessory parameter is involved. Performing the substitution z = R(t) ex- 
plicitly is useful. Substituting z = R(t) into (0) 'pulls it back' (cf. [TH]) to 



d 2 y ( R R r^^-npil^ abR2 n ^ 

^ + \-I + W^) [c - {a + b + m STt-w^R) y = - (3 - 3) 

To save space here, dR/dt, d 2 R/dt 2 are written as R, R. 

Lemma 3.4. The coefficient of the dy/dt term in the pulled-back hypergeo- 
metric equation ()3.3j) . which may be denoted W(t), equals the coefficient of the 
du/dt term in the Heun equation fiffi, i.e., •y/t + 6/(t — 1) + e/(t — d), if and 
only if R maps exponents to exponents. That is, the transformation at least 
partly 'works ' if and only if R maps exponents to exponents. 



PROOF. This follows by elementary, if tedious calculations. Suppose that 
R maps t = to to z = z = f R(t ) with multiplicity k, i.e., to leading order 
R(t) ~ Zq + C(t — t ) k ; if t and z are finite, that is. By direct computation, 
the leading behavior of W at t = to is the following. In the case when to is 
finite, W{t) ~ {l-k){t-to)- x if z ^ 0,1, oo; [1 - fc(l -c)](t-t ) _1 if Zo = 0; 
[1 — k(c — a — b)](t — to)' 1 if Zq = 1, and [1 — k(a + b)](t — to)' 1 if Zq = oo. 
In the case when t = oo, W(t) ~ (1 + k)t~ x if z ^ 0, 1, oo; [1 + k(l - c)}^ 1 
if zo = 0; [1 + k(c — a — b)}^ 1 if zo = 1, and [1 + k(a + b)]^ 1 if zo = oo. 

This may be restated as follows. At t = to, for finite to, the leading behavior 
of W is W(t) ~ (1— krf)(t— to)' 1 , where k is the multiplicity of t i— > z = f R(to) 
and 7] is the sum of the two characteristic exponents of the hypergeometric 
equation at z = z ; if the coefficient 1 — krj equals zero then W has no pole 
at t — t . Likewise, the leading behavior of W at t = oo is W{t) ~ (1 + krfjt' 1 , 
where k is the multiplicity of oo i— > z == R(oo) and rj is the sum of the 
two exponents at z = Zq', if the coefficient 1 + krj equals zero then W has a 
higher-order zero at t — oo. 
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By the definition of 'mapping exponents to exponents', it follows that the 
leading behavior of W at t = to, for all to finite, is of the form W(t) ~ 
(1 — i]')(t — to) -1 > and also at t = oo, is of the form W(t) ~ (1 + r)')t~ 1 , where 
in both cases r/ is the sum of the exponents of the Heun equation at t = t , 
iff R maps exponents to exponents. 

That is, the rational function W has leading behavior / yt~ 1 at t — 0, 5(t — 

at t — 1, e(t — cfj^ 1 at t = d, (1 + a + /3)t~ 1 = (7 + 5 + e)t _1 at t — 00, and is 

regular at all t other than 0, 1, d, 00, iff R maps exponents to exponents. □ 

The following two propositions characterize when the pulled-back hypergeo- 
metric equation ()3.3|) is, in fact, the Heun equation The first deals with 
Heun equations that are trivial in the sense of Definition 12. 1| and will be used 
in §0 The second will be applied to prove Theorem 13. II 

Proposition 3.5. A Heun equation of the form (H^J), which is trivial (i.e., 
a/3 = andq = 0), will be reduced to a hypergeometric equation of the form (RjJ 
by a specified rational substitution z = R(t) of its independent variable if and 
only if R maps exponents to exponents. 

PROOF. The 'if half is new, and requires proof. By Lemma l3.4[ the coef- 
ficients of dy/dt and du/dt agree iff R maps exponents to exponents, so it 
suffices to determine whether the coefficients of y and u agree. But by triv- 
iality, the coefficient of u in (^} is zero. Also, t = 00 has zero as one of its 
exponents, so all points t £ CP 1 have zero as an exponent. By the mapping 
of exponents to exponents, z = 00 must also have zero as an exponent, i.e., 
ab = 0. So the coefficient of y in ()3.3|) is also zero. □ 

Proposition 3.6. A Heun equation of the form (jffi), which has four singular 
points and is nontrivial (i.e., a[3 7^ or q 7^ 0), will be reduced to a hypergeo- 
metric equation of the form (ml by a specified rational substitution z = R(t) 
of its independent variable if and only if R maps exponents to exponents, and 
moreover, R is a polynomial, a/3 7^ ; and one of the following two conditions 
on the normalized accessory parameter p = q/a/3 is satisfied. 

(1) p equals one of 0,1, d. Call this point do, and the other two singular points 
&i and g?2- In this case, d must be a double zero of R or S, and each of 
c?i, g?2 must be a simple zero of R or S. 

(2) p does not equal any of 0,1, d. In this case, each of 0,1, d must be a 
simple zero of R or S , and p must be a triple zero of either R or S . 

In both cases, R and S must have no additional simple zeroes or zeroes of 
order greater than two. Also, if 1 — c (the nonzero exponent at z = 0) does not 
equal 1/2, then R must have no additional double zeroes; and if c — a — b (the 
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nonzero exponent at z = 1) does not equal 1/2, then S must have no additional 
double zeroes. 

Moreover, in both cases no additional double zero, if any, must be mapped 
by R to the point (out of z = 0, 1) to which p is mapped. (So additional double 
zeroes, if any, must all be zeroes of R, or all be zeroes of S .) 



PROOF. Like Proposition 13.51 this follows by comparing the pulled-back 
hypergeometric equation ()3.3|) to the Heun equation By Lemma EOl the 
coefficients of dy/dt and du/dt agree iff R maps exponents to exponents, so 
it suffices to characterize when the coefficients of y and u agree. 

The coefficient of y in (j3.3|) is to equal the coefficient of u in (J^J . It follows that 
ah = is possible iff a (3 = and q = 0, which is ruled out by nontriviality. So 
ah 7^ 0, and equality of the coefficients can hold iff 

= dR/dt dS/dt -(apt -q) fob = Co _Cj_ _C^_ 

~ R S t(t-l)(t-d) ~ t t-1 t-d' [ ' } 

where S = 1 — R, and at least two of Cq, Ci, G C are nonzero. 

Both Rr x dR/dt and S~ 1 dS/dt are sums of terms of the form n(t — A) -1 , 
where n is a nonzero integer and A is a zero or a pole of R or S. Poles are 
impossible, since A is a pole of R iff A is a pole of S, and there are no double 
poles on the right-hand side of ()3.4|) . So R must be a polynomial. 

By examining the definition of U in terms of R and S, one sees the following 
is true of any A G C: if R or S has a simple zero at t = A, then U will have a 
simple pole at t = A; if R or S has a double zero at t = A, then U will have 
an ordinary point (non-zero, non-pole) at t = A, and if R or S has a zero of 
order k > 2 at t = A, then U will have a zero of order — 2 at t = A. 

Most of what follows is devoted to proving the 'only if half of the proposition 
in the light of these facts, by examining the consequences of the equality (J3.4)) . 
In the final paragraph, the 'if half will be proved. 

There are exactly three ways in which the equality (|3.4|) can hold. 

(0) a/3 = 0, but due to nontriviality, q ^ 0. U has three simple poles on C, 
at t — 0, 1, d. It has no other poles, and no zeroes. So each of 0, 1, d must 
be a simple zero of either R or S; also, R and S can have no other simple 
zeroes, and no zeroes of order k > 2. Except for possible double zeroes, 
the zeroes of R and S are determined. The degree of R must equal the 
number of zeroes of R, and also equal the number of zeroes of S, counting 
multiplicity. But irrespective of how many double zeroes are assigned to 
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R or S, either R or S will have an odd number of zeroes, and the other 
an even number, counting multiplicity. So case cannot occur. 

(1) af3 7^ and a(3t — q is a nonzero multiple of t — d , where d = 0, 1, 
or d, so exactly one of C ,Ci,Cd is zero. U has two simple poles on C, 
at t = di,d2 (the two singular points other than do); it has no other 
poles, and no zeroes. So each of d\, c?2 must be a simple zero of either R 
or S; also, R and S can have no other simple zeroes, and no zeroes of 
order k > 2. Since by assumption (fj) has four singular points, each of 
0, 1, d must be a singular point, so the coefficient of dy/dt in (|3.3j) must 
have a pole at t — do, which implies that R or S must have a zero at do 
of the only remaining type: a double zero. 

(2) af3 7^ but a/3t—q is not a multiple of t, t—1, or t—d, so none of Co, C\, Cd 
is zero. U has three poles on C, and exactly one zero, at t = p = q/a/3, 
which is simple. So each of 0, 1, d must be a simple zero of either R or S, 
and must be a triple zero of either R or 5. Also, R and 5 can have 
no other simple zeroes, and no other zeroes of order k > 2. 

In cases 1,2, what remain to be determined are the (additional) double zeroes 
of R and S, if any. That is, it must be determined if any ordinary point of 
the Heun equation can be mapped to z = or z = 1 with double multiplicity. 
But by Proposition 13.31 R can map an ordinary point t = to to z = (resp. 
z — 1) in this way only if the exponents of z = (resp. z — 1) are 0, 1/2. 

Suppose this occurs. In case 1, if the exponents of t — p — do are denoted 0, 70, 
the exponents of R(p) will be 0, 7o/2, since t — p will be mapped with double 
multiplicity to z — R{p). So if R{to) = R(p) then 7 must equal 1, which, since 
q = aj3d , is ruled out by the assumption that each of 0, 1, d, including d , is 
a singular point. It follows that in case 1, R(t ) 7^ R{p)- A related argument 
applies in case 2. In case 2, the point p is an ordinary point of (5}), and a 
double critical point of the t 1— > z map. So as a singular point of (fj), R(p) 
must have exponents 0, 1/3. It follows that R(t ) = R{p) is impossible. 

The 'only if half of the proposition has now been proved; the 'if half remains. 
Just as ()3.4|) implies the stated conditions on R, so the stated conditions must 
be shown to imply (|3.4|) . But the conditions on R are equivalent to the left and 
right-hand sides having the same poles and zeroes, i.e., to their being the same 
up to a constant factor. To show the constant is unity, it is enough to consider 
the limit t — > 00. If degi? = n, then _R _1 d_R/dt ~ n/t and S~ 1 dS/dt ~ —n/t, 
so U, i.e., the left-hand side, has asymptotic behavior —n 2 /t 2 . This will be 
the same as that of the right-hand side if (af3)/(ab) = n 2 . But a = a/n and 
b = 13 /n follow from the assumption that R maps exponents to exponents. □ 



Finally, we can prove the main theorem, with the aid of the polynomial ma- 
nipulation facilities of the Macsyma computer algebra system. 
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Proof (of Theorem 13. 1|) By Proposition 13.61 the preimages of z = 0, 1 
under R must include t = 0,1, d, and in case 2, t = p = q/a(3. They may also 
include I (additional) double zeroes of R or of S, which will be denoted t = 
di, . . . , a\. Cases 1,2 of the theorem correspond to cases 1,2 of the proposition, 
and the subcases of the theorem correspond to distinct choices of I. 

Necessarily degR = |i? _1 (0)| = |i? _1 (l)|, where the inverse images are defined 
as multisets rather than sets, to take multiplicity into account. This places 
tight constraints on I, since each of 0, 1, d (and p, in case 2) may be assigned 
to either i? _1 (0) or i? _1 (l), but by the proposition, all of a 1; . . . , a/ must be 
assigned, twice, to one or the other. In case 1, one of 0,1, d (denoted do in 
the proposition) has multiplicity 2, and the other two (denoted dx,d 2 ) have 
multiplicity 1. It follows that < / < 2, with degR = I + 2. In case 2, each 
of 0, 1, d has multiplicity 1, and p has multiplicity 3. It follows that < / < 3, 
with deg R = I + 3. Subcases are as follows. 

(la) Case 1,1 — 0, degR = 2. Necessarily i? _1 (0), -R~ 1 (l) are {d , d }, {d\, d 2 }, 
or vice versa. Without loss of generality (w.l.o.g.), assume the latter, 
and assume d = 1. Then i? _1 (0) = {0,d} and i?" 1 (l) = {1,1}, i.e., 
S _1 (0) = {1, 1} and S _1 (l) = {0,d}. Since t — 1 is a double zero of S, 
S(t) = C{t - If for some C. But S(0) = 1, which implies C = 1, and 
S(d) = 1, which implies d = 2. So S(t) = (t - l) 2 and R(t) = t(2 - t). 
Since t = 0, d are both mapped singly to the singular point z = by R, 
their exponents must be those of z = 0, and hence must be identical. 

(lb) Case 1, I = 1, degR = 3. Necessarily _R _1 (0), i? _1 (l) are {d ,do,di}, 
{d%, a\,ax\, or vice versa. W.l.o.g., assume the former, and also assume 
do, di, d2 equal 1, d, 0, respectively. Then R~ l {0) = {1, 1, d} and i? _1 (l) = 
{0, ax j a i}- It follows that R(t) — (t — 1) 2 (1 —t/d), where d is determined 
by the condition that the critical point of R other than t = 1 (i.e., t = 
ai) be mapped to 1. Solving dR/dt = yields a\ = (2d + l)/3, and 
substitution into R(ai) — 1 = yields d = 4 or —1/2. But the latter 
is ruled out: it would imply a\ = 0, which is impossible. So d = 4 and 
a\ = 3. Since t = l,d are mapped to the singular point z = 0, doubly and 
singly respectively, the exponents of t = 1 must be twice those of z = 0, 
and the exponents of £ = d must be the same as those of z = 0. 

(lc) Case 1,1 — 2, degR = 4. Necessarily i? _1 (0) and i? _1 (l) are {d , d , di, d 2 } 
and {ax, ax, a 2 , a 2 }, or vice versa. W.l.o.g., assume the latter, and as- 
sume do = 1. Then i? _1 (0) = {ax, ax, a 2 , a 2 } and i? _1 (l) = {0,1, l,d}, 
i.e., S'~ 1 (0) = {0, 1,1, d} and 5 ,_1 (1) = {ai,ai,a 2 ,a 2 }. So S(t) equals 
t(t — l) 2 (t — d), where d is determined by the condition that S must have 
two critical points other than t — 1, i.e., t = a\,a 2 , which are mapped 
by S to the same critical value (in fact, to z = 1). Computation yields 
dS/dt = A(t — 1) [4t 2 — (3d + 2)t + d], so ai,a 2 must be the roots of 
4t 2 — (3d + 2)t + d. If the corresponding critical values are Aw\, Aw 2 , then 
Wi,w 2 are the roots of the polynomial in w obtained by eliminating t 
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between w — S(t)/A and 4t 2 — (3d + 2)t + d. Its discriminant turns out 
to be proportional to (d — 2) 2 (9d 2 — Ad + 4) 3 , so the criterion for equal 
values is that d = 2 or 9c? 2 — Ad + A = 0. But the latter can be ruled out, 
since by examination it would result in a±,a 2 being equal. So d = 2; 

ai,a 2 = 1 ± v^/2; and = At(t - lf(t - 2) with A = -4, so that 

r 1 1 2 

S(ai) = 1. Hence = 4 t(2 — t) — | . Since t — 0, d are mapped 
simply to z — 1 and t = 1 is mapped doubly, the exponents of t — 0, d 
must be the same, and double those of t — 1. 

(2a) Case 2, / = 0, deg.R = 3. Necessarily iT^O), are {0, 1, d}, 

or vice versa. W.l.o.g., assume the former. Then R(t) = A(t — p) 3 for 
some A. Since t = 0, 1, d are to be mapped singly to 1, they must be 
the vertices of an equilateral triangle, with mean p, so A = — 1/p 3 and 
R — (1 —t/p) 3 . W.l.o.g., take d — \ + i^, so p = \ + i^. The exponents 
at t — 0, 1, d must be equal, since they all equal the exponents at z — 1. 

(2b) Case 2, / = 1, degi? = 4. Assume w.l.o.g. that i? _1 (0) = 6?} 
and = {0,l,ai,ai}, i.e., 5-^0) = {0,l,ai,ai} and ^^l) = 

{p,p,p, d}. It follows that R(t) = (1 — t/d)(l — t/p) 3 , but to determine g? 
and p, it is best to focus on S. Necessarily S(t) = At(t — l)(t — 6*i) 2 , and 
p can be a triple zero of R iff it is a double critical point of S as well as R. 
The condition that S have a double critical point determines a\. dS/dt = 
A(t - a ± ) [At 2 - (3 + 2ai)t + ai], so the polynomial 4t 2 - (3 + 2a x )* + a x 
must have a double root. Its discriminant is 4a 2 — 4ai + 9, which will 
equal zero iff a± — | ± i\/2. The corresponding value of the double root, 

i.e., the mandatory value of p, is | ± i^. The requirement that S map 
p to 1 implies A = l/p(p — l)(p — ai) 2 . d is determined as the root 
of R = 1 — S other than p; some computation yields \ ± i^j^. W.l.o.g. 
the '±' in the expressions for p and d can be replaced by '+'. Since t — p 
is an ordinary point and R maps t = p triply to z = 0, z = must have 
exponents 0, 1/3. Since R maps t — d simply to z — 0, t — d must also 
have exponents 0, 1/3. Similarly, since R maps the ordinary point t — a\ 
doubly to z — 1, z — 1 must have exponents 0, 1/2; so t = and t — 1, 
which are mapped simply to z = 1, must also. 
(2c) Case 2, I = 2, degi? = 5. Assume w.l.o.g. that i2 _1 (0) = {p,p,p,0, 1} 
and -R _1 (l) = {d, ai, ai, a 2 , a 2 }. Then i?(t) = At(t — 1)(£ — p) 3 , where 
p is determined by i? having two critical points other than t = p, i.e., 
t = ai,a 2 , which are mapped to the same critical value (i.e., to z — 1). 
dR/dt = A(t — p) 2 [5t 2 — (2p + A)t + p], so a\, a 2 must be the two roots 
of 5t 2 — (2p + 4)t + p. If the corresponding critical values are Aw±, Aw 2 , 
then w±,w 2 are the roots of the polynomial in w obtained by eliminating t 
between w — R(t)/A and 5t 2 — (2p + 4)t + p. Its discriminant turns out 
to be proportional to (p 2 — p + 4) 3 (27p 2 — 27p + 8) 2 , so the criterion for 
equal values is that p 2 — p + 4 = or 27p 2 — 27p + 8 = 0. But the former 
can be ruled out, since by examination it would result in a±,a 2 being 
equal. The latter is true iff p — \ ± i^jp- W.l.o.g., the plus sign may be 
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used. This yields ai, a 2 = ~ ± 2 -^p + i^p- From the condition i?(a*) = 1, 
it follows that A = — i 202 ^ 1 ^ . d is determined as the root of - 1 
other than ai,a 2 ; computation yields d = | + ■ Since t = p is an 
ordinary point mapped triply to z = 0, z = must have exponents 0, 1/3. 
Similarly, since R maps the ordinary points t = di to z = 1, z = 1 must 
have exponents 0, 1/2, so t = d, which is mapped singly to it, must also. 
(2d) Case 2, 1 = 3, deg R = 6. Necessarily iT^O) and i? _1 (l) are {p,p,p, 0, 1, d} 
and {ai, ai, 02, a 2 , 03, a 3 }, or vice versa. W.l.o.g., assume the latter. Then 
R(t) = A(t - ai) 2 {t - a 2 f{t - a 3 ) 2 and S(t) = Bt(t - l)(t - d)(t - p) 3 . 
Since t = p is a triple zero of S, R(t) ~ 1 — C(t — p) 3 for some nonzero C. 
So \J R{t) 1 defined to equal +1 at t — p, will have a similar Taylor series: 

VA(t — a\)(t — a 2 )(t — a 3 ) ~ 1 — C(t — p) 3 /2. This is possible only if 
ai,a 2 ,a 3 are the vertices of an equilateral triangle, and p is their mean. 
It follows that the roots of S other than t = p, i.e., t = 0,1, d, are 
also the vertices of an equilateral triangle centered on p. W.l.o.g., choose 
d — I + i^ and p = | + i^p With a bit of algebra, R can be rewritten 
in the form given in the theorem. Since t = p is an ordinary point and 
R maps it triply to z = 0, z = must have exponents 0, 1/3. Since R maps 
t = 0, 1, d simply to z = 0, t = 0, 1, d must also have exponents 0, 1/3. 

The theorem is proved. □ 



Now that the main theorem is proved, we can proceed to derive explicit Heun- 
to-hypergeometric reduction formulae. Theorem 13 . 71 gives a necessary condition 
for the existence of a reduction, and Theorem 13 . Sl presents a representative list. 

Theorem 3.7. Suppose the Heun equation (^J) has four singular points and 
is nontrivial (a/3 7^ or q 7^ 0). Then its local solution Hl(d,q;a,/3,j,d;t) 
can be reduced to a hypergeometric function 2 Fi(a, b; c; z) by a formula of the 
type Hl(t) = 2 Fi(R(t)), with R a rational function, only if its parameters 
d, q; a, [3 satisfy q = af3p, with (d,p) equal to one of the following 23 pairs. If 
a reduction of this type exists, R(t) will be a polynomial of the stated degree. 



(la) [deg R = 2 or 4]. 

(lb) [degR = 3]. 

(2a) [deg R = 3 or 6]. 

(2b) [degR = 4]. 

(2c) [degi2 = 5]. 
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PROOF. The five subcases are taken from Theorem 13.11 with (lc),(2d) re- 
spectively subsumed into (la), (2a). Theorem 13.11 supplies only a single pair 
(D,po) for each subcase. By the remarks preceding that theorem, d may be 
any value on the cross-ratio orbit of D. There are up to 6 distinct possibilities, 
which were given in ([2.7)1 . The value of p associated to d is the preimage of po 
under the corresponding affirie map 1 1— > Ai(t), given in ()2.8|) . □ 

Theorem 3.8. Suppose a Heun equation has four singular points and is non- 
trivial (a/3 or q ^ 0). Then the only reductions of its local Heun func- 
tion HI to 2 i*i that can be performed by a rational transformation of the in- 
dependent variable involve polynomial transformations of degrees 2, 3, 4, 5, 
and 6. There are seven distinct types, each of which can exist only if d lies on 
an appropriate cross-ratio orbit. The following list includes a representative 
reduction of each type. The ones with real d {and deg R = 2, 3, 4) include 



HI (2, a/3; a, /3, 7, a + (3 - 27 + 1; t) 

= 2 Fx{% f; 7; *(2-t)), 
HI (4, a(3; a, (3, §, %2±&; t) 

= 2 Fi (f, f; |; i-(t-i) 2 (i-t/4) 



JH(2, a/9; a, /3, 2±|±2, 2±5; t 



a /?. a+/3+2. 



1 1 41 4 ' 



t(2 -*) 



(3.5a) 
(3.5b) 
(3.5c) 



and the ones with non-real d (and deg R = 3,4, 5, 6) include 



Hl(±±i&,a(3(i±if);a, (3, 



Q+/3+1 a+8+1. 



^1 



3' 3' 3 ' 



1 - 



*/(§±if) 
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a 1 
4' 6 



a. 1. 

4' 2' 



l-t/(i±i¥) l-*/(i±i^) 



(2 ^ 90 ' ^"(e < ^)(2~^^18^' < ^''6 ^'3' 3' ^ 



90 

pin 1 - 

21 I 5 ' 6 5 ' 3' 



a. 2. 



(Ti 



2025V15 
64 



*(*-!) 



18 

(I ± , a(l - a) (| ± i^); a, 1 - a, §, §, t 



,Fi I - - - -• 

- 1 I 6 ' 6 fi ' 



a. 2. 

6 ' 3' 



1-4 



l-t/(i±if 



(3.6a) 



(3.6b) 



(3.6c) 



(3.6d) 



In the preceding reductions, a, (3, 7 are free parameters. Each of these equalities 
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holds in a neighborhood of t = whenever the two sides are defined, e.g., 
whenever the fifth argument of HI and the third argument 0/2 -Pi are not equal 
to a nonpositive integer. 

Remark 3.8.1. The equalities of Theorem 13.81 hold even if the Heun equation 
has fewer than four singular points, or is trivial; but in either of those cases, 
additional reductions are possible. For the trivial case, see §0 

Remark 3.8.2. The special harmonic reduction (|3.5cj) is composite: it can be 
obtained from the case 7 = (a + (3 + 2) /4 of the reduction (|3.5aj) by applying 
Gauss's quadratic hypergeometric transformation (21 §3.1] 

2 F 1 (a, b; (o+6+l)/2; z) = 2 F l (a/2, b/2; (a+6+l)/2; l-A(z-\) 2 ) (3.7) 

to the right-hand side. The special equianharmonic reduction (|3.6d|) can be 
obtained in the same way from the case f3 = 1 — a of the reduction ()3.6a|) . 

One might think that by applying ()3.7|) to the right-hand sides of the remaining 
reductions in (J3.5a|) - (|3.6d|) . additional composite reduction formulae could be 
generated. However, there are only a few cases in which it can be applied; and 
when it can, it imposes conditions on the parameters of HI which require that 
the corresponding Heun equation have fewer than four singular points. 

PROOF. HI and 2 -Pi are the local solutions of their respective equations 
which belong to the exponent zero at t — (resp. z — 0), and are regular and 
normalized to unity there. So the theorem follows readily from Theorem 13.11 
l)3.5aj) - (|3.5c|) come from subcases la-lc, and (|3.6ajl — (|3.6djl from subcases 2a- 
2d. In each subcase, the Gauss parameters (a, b; c) of 2 -Pi are computed by first 
calculating the exponents at z — 0, 1, 00, in the way explained in Remark 13.1.31 
In some subcases, the polynomial map supplied in Theorem 13.11 must be chosen 
to be S = 1 — R rather than R, due to the need to map t = to z = rather 
than to z = 1, so that the transformation will reduce HI to 2 -Pi, and not to 
another local solution of the hypergeometric equation. □ 

The list of Heun-to-hypergeometric reductions given in Theorem 13.81 is repre- 
sentative rather than exhaustive. For each subcase of Theorem 13 .1\ there is 
one reduction for each allowed value of d. Each reduction on the above list 
came from choosing d = D, but any other d on the cross-ratio orbit of D may 
be chosen. The orbit is defined by AOld being one of the triangles (at most 
six) similar to A01D, i.e., by A01Z) being obtained from AOld by an affine 
transformation A\ G A(Sj). So for any subcase of Theorem 13.11 and choice 
of d, the appropriate polynomial map will be z = A 2 (Ri(Ai(t))) , where Ai is 
constrained to map AOlc? to AO 1£) and is listed in (|2.8|) . Ri is the polynomial 
map given in the subcase, and A 2 G A(§), i.e., A 2 {z) = z or 1 — z, is chosen 
so that t = is mapped to z = rather than to z — 1. 



22 



As an example, consider the harmonic subcase la of Theorem 13.11 in which 
D = 2, the cross-ratio orbit of D is { — 1, f,2}, and the polynomial map is 
R\{t) — t(2 — t). Choosing d = D yields the reduction (|3.5a|) . Choosing d = 
1 — D = — 1 yields an alternative reduction of HI to 2 Fi, namely 

Hl(-1, 0; a, (3, 7 , (a + (3 - 7 + l)/2; t) = 2 F X (a/2, (3/2; ( 7 + l)/2; t 2 ) , 

(3.8) 

in which Ai(t) — 1 — t according to (|2.8|l . and ^(z) = 1 — 2. 

It is not difficult to check that in all, exactly 28 Heun-to-hypergeometric re- 
ductions can be derived from Theorem 13.11 They exhibit the 23 values of the 
pair (d,q/af3) listed in Theorem 13.71 Of the 28, eleven were given in Theo- 
rem EH1 and ()3.8|) is a twelfth. With the exception of the two reductions with 
(d,q/afl) = (—1,0), one of which is (|3.8|) . the 28 split into pairs, each pair 
being related by the identity (|2.1U|) . which takes d to 1/d. 



4 A Generalization 



In applied mathematics, it is seldom the case that the four singular points 
of an equation of Heun type are located at 0, 1, d, oo. But the main theorem, 
Theorem 13.11 may readily be generalized. Consider the situation when three 
of the four have zero as a characteristic exponent, since this may always be 
arranged by applying an F-homotopy. There are two cases of interest: either 
the singular points include oo and each of the finite singular points has zero as 
a characteristic exponent; or the location of the singular points is unrestricted. 
The latter includes the former. They have respective P-symbols 





di 


da 


d 3 


oo 




d\ d 2 d 3 


di 


p < 











a ; s 







a ; s 




. !"7 


1-5 


1-e 






1-7 1-6 1-e 


P 



(4-1) 

In the nomenclature of Ref. |22j . they are canonical cases of the natural 
general-form and general-form Heun equations. They are transformed to the 
Heun equation (S)) by the affine and Mobius transformations 



s — di 
d 2 — di 



(s - d 1 )(d 2 - d 4 ) 
(d 2 - di)(s - d^ 



(4.2) 
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respectively. Each of the P-symbols (j4.1j) is accompanied by an accessory 
parameter. The equation specified by (|4.1b ) can be written as 




a(3s -q' 

u = 0, 



ds 2 \s — di s — d 2 s — d 3 J ds (s — di)(s — d 2 )(s — d 3 ) 

(4.3) 

where q' is the accessory parameter j22|. The equation specified by (|4.1b ) with 

<J 4 ^ cx) can be written as 

d 2 u ( 7 5 e 1 — a — d\ du 

3 



a/3 

+ 



i=l 




11(^4 - di) / (s-d 4 )-q" 

u = 0, (4.4) 



(s — di)(s — d 2 )(s — dz)(s — d 4 
where q" is the accessory parameter 



The impediment to the generalization of Theorem 13. II to these two equations 
is the specification of the cases that should be excluded due to their being 
'trivial', or having fewer than four singular points. The excluded cases should 
really be specified not in terms of the ad hoc parameters q' and q", but rather 
in an invariant way, in terms of an accessory parameter defined so as to be 
invariant under affine or Mobius transformations, respectively. Schafke |2Sj has 
defined new accessory parameters of second-order Fuchsian equations on CP 1 
that are invariant under affine transformations, but no extension to general 
Mobius transformations seems to have been developed. 

In the absence of an invariantly defined accessory parameter, an ad hoc ap- 
proach will be followed. It is clear that (J4.3|) is trivial, i.e., can be transformed 
to a trivial Heun equation by an affine transformation, iff a(3 — 0, q' — 0. Also, 
it will have fewer than four singular points if 7 = 0, q' = 0; or 5 = 0, q' = a/3; 
or e = 0, q' = af3d. Likewise, it is fairly clear that (j4.4j) will be trivial, i.e., 
can be transformed to a trivial Heun equation by a Mobius transformation, 
iff a/3 = 0, q" = 0. The conditions on the parameters for there to be a full set 
of singular points are, however, more complicated. 

The first generalization of Theorem 13.11 is Corollary 14. 1[ which follows from 
Theorem I3.1l bv applying the affine transformation (|4.2b ). It mentions a poly- 
nomial transformation, which is the composition of the s 1— > t affine transfor- 
mation with the 1 1 — ► z polynomial map of Theorem 13.11 To avoid repetition, 
Corollary 14.11 simply cites Theorem 13.11 for the necessary and sufficient condi- 
tions on the exponent parameters and the accessory parameter. 

Corollary 4.1. A natural general-form Heun equation of the canonical type 
(14. 3|) . which has four singular points and is nontrivial {i.e., a/3 7^ or q' 7^ 0), 
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can be reduced to a hypergeometric equation of the form ^ by a rational 
substitution z = R{s) iff a(3 ^ 0, R is a polynomial, and the Heun equation 
satisfies the following conditions. 

(i) Adid 2 dz must be similar to A01D, with D = 2 or ~ + or D = 4, 
i_l_j5^|^ or i _|_- ii\/i5 _ j 1 /^ ^ ^ mwsi either be a degenerate triangle consisting 
of three equally spaced collinear points {the harmonic case), or be an equilateral 
triangle [the equianharmonic case), or be similar to one of three other specified 
triangles, of which one is degenerate and two are isosceles, (ii) The exponent 
parameters 7, 5, e must satisfy conditions that follow from the corresponding 
subcases of Theorem 13. 1L (iii) The parameter q' must take a value that can be 
computed uniquely from the parameters 7, 5, e and the choice of subcase. 

Example 4.1.1. In the harmonic case, the two endpoints of the degenerate 
triangle of singular points Adid 2 d% must have equal exponent parameters, 
and q' must equal a/3 times the intermediate point. In this case, R will typ- 
ically be a quadratic polynomial. There are two possibilities: R will map the 
two endpoints to z = and the intermediate point to z — 1, or vice versa. 
If the characteristic exponents of the intermediate point are twice those of the 
endpoints, then R may be quartic instead: the composition of either possible 
quadratic polynomial with a subsequent z 1— > A(z — |) 2 or z 1— > Az(l — z) map. 

Example 4.1.2. In the equianharmonic case, all three exponent parameters 
7, 5, e must be equal, and the accessory parameter q' must equal a(3 times the 
mean of di, d 2 , d 3 . In this case, R will typically be a cubic polynomial. There 
are two possibilities: R will map d±, d 2 , d^ to z = and their mean to z = 1, or 
vice versa. If the exponent parameters 7, 8, e equal 2/3, then R may be sextic 
instead: the composition of either possible cubic polynomial with a subsequent 
z 1— > A(z — |) 2 or z t— > Az(l — z) map. 

The further generalization of Theorem l3.1l is Corollary 14. 21 which follows from 
Theorem 13.11 by applying the Mobius transformation (|4.2b ). It mentions a 
rational substitution, which is the composition of the s 1— > t Mobius transfor- 
mation with the 1 1— > z polynomial map of Theorem 13.11 

Corollary 4.2. A general-form Heun equation of the canonical type (|4.4jl . 
which has four singular points and is nontrivial (i.e., a/3 7^ or q" 7^ 0), can be 
reduced to a hypergeometric equation of the form ^ by a rational substitution 
z = R{s) iff a(3 7^ 0, and the Heun equation satisfies the following conditions. 

(i) The cross-ratio orbit of {di, d 2 , c? 3 , d^} must be that of {0, 1,D, oo} ; where 
D is one of the five values enumerated above. That is, it must be the harmonic 
orbit, the equianharmonic orbit, or one of three specified generic orbits, one 
real and two non-real, (ii) The exponent parameters 7, 8, e must satisfy con- 
ditions that follow from the corresponding subcases of Theorem 13.11 (iii) The 
parameter q" must take a value that can be computed uniquely from the pa- 
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rameters 7, 5, e and the choice of subcase. 

Example 4.2.1. Suppose di,d 2 , ^3,^4 form a harmonic quadruple, i.e., can be 
mapped by a Mobius transformation to the vertices of a square in C. More- 
over, suppose two of d±,d 2 , d% have the same characteristic exponents, and are 
mapped to diagonally opposite vertices of the square. That is, of the three 
parameters 7, 5, e, the two corresponding to a diagonally opposite pair must 
be equal. Then provided q" takes a value that can be computed from the other 
parameters, a substitution R will exist. It will typically be a degree-2 rational 
function, the only critical points of which are the third singular point (out of 
di,d2,d 3 ) and d^. Either R will map the two distinguished singular points to 
z = 1 and the third singular point to z = 0, or vice versa; and d^ to z = 00. 
In the special case when the exponents of the third point are twice those of the 
two distinguished points, it is possible for R to be a degree-4 rational function. 

Example 4.2.2. Suppose d±, d 2 , d 3 , d± form an equianharmonic quadruple, i.e., 
can be mapped by a Mobius transformation to the vertices of a regular tetra- 
hedron in CP 1 . Moreover, suppose di,d 2 ,d 3 have the same characteristic ex- 
ponents, i.e., 7 = 5 = e. Then provided q" takes a value uniquely determined 
by the other parameters, a substitution R will exist. Typically, R will be a 
degree-3 rational function, the only critical points of which are the mean of 
di,d 2 ,d 3 with respect to d±, and dj. Either R will map di,d 2 ,d 3 to z = 1 
and the mean of di,d 2 ,d% with respect to dj to z = 0, or vice versa; and d^ 
to z = 00. In the special case when the exponents of each of di,d 2 ,d 3 equal 
0, 1/3, it is possible for R to be a degree-6 rational function. 

Remark 4.2.3. In Example 14.2.21 the concept of the mean of three points 
in CP 1 with respect to a distinct fourth point was used. A projectively invari- 
ant definition is the following. If T is a Mobius transformation that takes d& 
(7^ di,d 2 , g? 3 ) to the point at infinity, the mean of di,d 2 , d 3 with respect to d^ 
is the point that would be mapped to the mean of Td±,Td 2 , Td% by T. 



5 The Clarkson— Olver Transformation 

The reduction discovered by Clarkson and Olver |6 a , which stimulated these 
investigations, turns out to be a special case of the equianharmonic Heun- 
to-hypergeometric reduction of §H3 Their reduction was originally given in a 
rather complicated form, which we shall simplify. 

Recall that the Weierstrass function p(u) = p(u; g 2 , #3) with invariants g 2 , g 3 G 
C, which cannot both equal zero, has a double pole at u = and satisfies 

p 2 = 4p 3 -g 2 p-g3 (5.1) 
= 4(p - ei)(p - e 2 )(p - e 3 ). 
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Here ei, e 2} e3, the zeroes of the defining cubic polynomial, are the finite critical 
values of p, the sum of which is zero; they are required to be distinct, p is 
doubly periodic on C, with periods denoted 2u, 2a/. So it can be viewed as a 
function on the torus T = C/£, where £ = 2u;Z © 2o/Z is the period lattice. 
It turns out that the half-lattice {0, u, a/, a; + a/} + £ comprises the critical 
points of p. The map p : T — > CP 1 is a double branched cover of the Riemann 
sphere, but T is uniquely coordinatized by the pair (p, p'). 

The modular discriminant A = f g\ — Hlg\ ^ is familiar from elliptic function 
theory. If g%,gz G K and A > (the so-called real rectangular case, which 
predominates in applications), u,u' can be chosen to be real and imaginary, 
respectively. If A < (the less familiar real rhombic case) , they can be chosen 
to be complex conjugates, so that the third critical point lu 2 = u> + a/ is real. 

Clarkson and Olver considered the Weierstrass-form Lame equation 

d 2 ^ 



du 2 



+ l)p(u) + B],p = ^ (5.2) 



which is a Fuchsian equation on T with exactly one singular point (at (p, p') = 
(oo, oo)) and a single accessory parameter, B. [We have altered their exponent 
parameter —36a to £{£ + 1), to agree with the literature, and have added the 
accessory parameter.] In particular, they considered the case g 2 = 0, g% ^ 0, 
B = 0. They mapped u G T to z G CP 1 by the formal substitution 

- 1 [ {1 - Z)1/S dT k * 

and showed that the Lame equation is reduced to 

^-^d? + (2-?)^ + 36 ^ = (5 ' 4) 

This is a hypergeometric equation with (a, b; c) = (—£/6, {£ + l)/6; 1/2). It has 
exponents 0, 1/2 at z = 0; 0, 1/3 at z = 1; and —£/6, {£ + l)/6 at z — oo. 

In elliptic function theory the case g 2 = 0, g% ^ is called equianharmonic, 
since the corresponding critical values e±, are the vertices of an equilateral 
triangle in C. If, for example, g% G M, then A < 0; and by convention, e±, e2, 
correspond to u,u 2 ,uj', respectively. e\ and are complex conjugates, and 
e 2 is real. The triangle AOo^a/ is also equilateral PQ § 18.13]. 

So, what Clarkson and Olver considered was the equianharmonic Lame equa- 
tion, the natural domain of definition of which is a torus T (i.e., a complex 
elliptic curve) with special symmetries. For the Lame equation ()5.2j) to be 
viewed as a Heun equation on CP 1 , it must be transformed by s = p(u) to its 
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algebraic form [13]. The algebraic form is 

f± + ( 1/2 + J/2_ + _l/2_\ dip + [-i(j + l)/A]s-B/A = Q 
ds 2 \s - ei s - e 2 s - e 3 / ds (s - e{)(s - e 2 )(s — e 3 ) 

(5.5) 

This is a special case of (|4.3|) , the canonical version of the natural general- form 
Heun equation, with distinct finite singular points dx,d 2 ,d 3 = ei,e 2 ,e 3 . Also, 
a,f3 = -1/2, (£ + l)/2, 7 = 8 = e = 1/2, and g' = 5/4. It has characteristic 
exponents 0, 1/2 at s = ex, e 2 , e 3 , and —£/2, (£ + l)/2 at s — oo. 

Applying Corollary 14. II to ()5.5|) yields the following. 

Theorem 5.1. The algebraic-form Lame equation ()5.5jl . m i/ie equianhar- 
monic case g 2 = 0, g% ^ 0, can be reduced when £(£ + 1) ^ to a hypergeo- 
metric equation of the form ^ by a rational transformation z = R(s) iff the 
accessory parameter B equals zero. If this is the case, R will necessarily be a 
cubic polynomial; 

z = 4s 3 /g 3 , z = l- 4s 3 /g 3 (5.6) 
will both work, and they are the only possibilities. 

PROOF. If £{£ + 1) ^ 0, the Heun equation (|5.5J) is nontrivial in the sense 
of Definition 12.11 with four singular points; by (|5.1jh the are the cube roots 
of #3/4, and are the vertices of an equilateral triangle. Since 7 = 5 = e, the 
equianharmonic case of Corollary 14.11 applies, and no other. 

The mean of ei, e 2 , e 3 is zero. So the polynomial 4s 3 /g 3 is the cubic polynomial 
that maps each singular point to 1, and their mean to zero; 1 — 4s 3 /g 3 does the 
reverse. These are the only possibilities for the map s h z, since the sextic 
polynomials mentioned in the equianharmonic case of Corollary 14.11 can be 
employed only if 7, 5, e equal 2/3, which is not the case here. □ 

Remark. Corollary 14.11 (equianharmonic case) also applies to the equianhar- 
monic algebraic-form Lame equation with £(£+ 1) = 0, B ^ 0, and guarantees 
it cannot be transformed to the hypergeometric equation by any rational sub- 
stitution; since in the sense used above, this too is a nontrivial Heun equation. 

Corollary 5.2. The Weierstrass-form Lame equation ()5.2j) . in the equianhar- 
monic case g 2 = 0, g 3 7^ ; can be reduced when £(£+1) 7^ to a hypergeometric 
equation of the form ^ by a substitution of the form z = R (p(u)), where R is 
rational, iff the accessory parameter B equals zero. In this case, 

z = Ap(u) 3 /g 3 , z = l- 4p(u) 3 /g 3 (5.7) 

will both work, and they are the only such substitutions. 
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Applying the substitution z = 1 — 4p(u) 3 /g 3 to the Lame equation ()5.2|) re- 
duces it to the hypergeometric equation ([5.4)1 . as is readily verified. (The other 
substitution z = 4p(u) 3 /g 3 yields a closely related hypergeometric equation, 
with the singular points z = 0, 1 interchanged.) The Clarkson-Olver substi- 
tution formula ()5.3|) contains a multivalued elliptic integral, but it may be 
inverted with the aid of (|5.1|) to yield z — 1 — Ap(u) 3 '/ ' g$. So their transforma- 
tion fits into the framework of Corollary 15.21 

The most noteworthy feature of the Clarkson-Olver transformation is that it 
can be performed irrespective of the choice of exponent parameter t. Only the 
accessory parameter B is restricted. As they remark, when i = 1, 1/2, 1/4, 
or 1/10, it is a classical result of Schwarz that all solutions of the hypergeo- 
metric equation ()5.4|) are necessarily algebraic fSJ §10.3]. This implies that 
if B = 0, the same is true of all solutions of the algebraic Lame equation ([5.5)1 : 
which had previously been proved by Baldassarri [5], using rather different 
techniques. But irrespective of the choice of £, the solutions of the B = 
Lame equation reduce to solutions of the hypergeometric equation. This is 
quite unlike the other known classes of exact solutions of the Lame equation, 
which restrict £ to take values in a discrete set [20, § 2.8.4]. But it is typical of 
hypergeometric reductions of the Heun equation. As the theorems of §01 make 
clear, in general it is possible to alter characteristic exponents continuously, 
without affecting the existence of a reduction to the hypergeometric equation. 

It should be mentioned that the harmonic as well as the equianharmonic case 
of Corollarv l4.1l can be applied to the algebraic-form Lame equation. One of the 
resulting quadratic transformations was recently rediscovered by Ivanov |16j . 
in a heavily disguised form. The case of quadratic rather than cubic changes 
of the independent variable will be considered elsewhere. 



6 The Seemingly Trivial Case a(3 = 0, q = 



If the Heun equation (^} is trivial in the sense of Definition 12.11 it may be 
solved by quadratures. A basis of solutions is 



ut(t) 



u 2 (t) 



cxp 




dv. (6.1) 



In the trivial limit, the local Heun function Hl(d, q; a, j3, 7, 5; t) degenerates to 
the former, and the solution belonging to the exponent 1 — 7 at t — 0, de- 
noted Hl(d, q; a, /3, 7, 5; t) here, to the latter. In applications, explicit solutions, 
if any, are what matter. It is nonetheless interesting to examine under what 
circumstances a trivial Heun equation can be reduced to a hypergeometric 
equation. This question was first considered by Kuiken [TSj . 
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The canonical polynomial substitutions of §HI]give rise to many nonpolynomial 
rational reductions of trivial Heun equations to hypergeometric equations, by 
composing with certain Mobius transformations. To understand why, recall 
that Theorem 13.11 characterized, up to affine automorphisms of the two equa- 
tions, the polynomial substitutions that can reduce a nontrivial Heun equation 
to a hypergeometric equation. If t i— > Ri(t) denotes a canonical polynomial 
substitution, the full set of polynomial substitutions derived from it comprises 
all t i— > A 2 (Ri(Ai(t))), where A\ G A(fi) is an affine automorphism of the 
Heun equation, which maps {0, 1, d} onto {0, 1, D}, and A 2 G A(t)) is an affine 
automorphism of the hypergeometric equation, which maps {0, 1} onto {0, 1}. 
(The only two possibilities are A 2 (z) = z and A 2 (z) — 1 — z.) 

In the context of nontrivial Heun equations, Mobius automorphisms that are 
not affine could not be employed; essentially because, as discussed in § 12.31 
moving the point at infinity would require a compensating F-homotopy. But 
in the trivial case no such issue arises: by Proposition 13.51 the Heun equation 
is reduced to a hypergeometric equation by a rational substitution of its inde- 
pendent variable, z = R(t), iff the substitution maps exponents to exponents. 
And Mobius transformations that are not affine certainly preserve exponents. 

Theorem 6.1. A Heun equation of the form flfi]), which has four singular 
points and is trivial (i.e., a(3 = and q = 0), can be reduced to a hyper- 
geometric equation of the form ((Q) by any rational substitution of the form 
z = M 2 (R\(Mi(t))), where z = Ri(t) is a polynomial that maps {0,1, D} 
to {0, 1} ; listed (along with D) in one of the seven subcases of Theorem 13.11 
and where Mi G Ji4(S}), and M\ G -M((}). That is, Mi maps {0, 1, d, oo} onto 
{0, 1, D, oo}, and M 2 maps {0, 1, oo} onto {0, 1, oo}. The necessary conditions 
on characteristic exponents stated in Theorem 13.11 must be satisfied, the condi- 
tions on exponents at specified values oft being taken to refer to the exponents 
at the preimages of these points under Mi . 

Remark 6.1.1. As in the derivation of the Hl(t) = 2 Fi(R(t)) reduction for- 
mulae listed in Theorem 13.81 the Gauss parameters (a, b; c) of the resulting 
hypergeometric equation can be computed by first calculating the exponents 
at z — R(t) = 0, 1, oo, using the mapping of exponents to exponents. 

The following example shows how such nonpolynomial rational substitutions 
are constructed. In the harmonic subcase la of Theorem 13.11 D = 2 and the 
polynomial tranformation is t t— > z = Ri(t) = t(2 — t); the necessary condition 
on exponents is that t — 0, d have identical exponents. Consider d = — 1, 
which is on the cross-ratio orbit of D. Mi(t) = (t — l)/t can be chosen; also, 
let M 2 (z) = 1/z. Then the composition 

z = R(t) = M 2 (i?i(Mi(i))) = t 2 /(t 2 - 1) (6.2) 

maps t = to z = and t = oo to z = 1 (both with double multiplicity), 
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and t = l,d to z = oo. This substitution may be applied to any trivial Heun 
equation with d = — 1 and identical exponents at t = 1, d, i.e., with 5 = e. 

In this example, Mi,M% were selected with foresight, to ensure that R maps 
t — to z — 0. This_makes it possible to regard the substitution as a reduction 
of HI to 2 or of HI to 2 -Pi- By calculation of exponents, the reduction is 

HI {-I, 0;0, /?, 7 ,(l + /? + 7)/2;t) (6.3) 
= (-l)^ 1 )/ 2 2^ (0, (1 - (5 + 7 )/2; (1 + 7 )A t 2 /(* 2 - 1)) • 

The normalization factor (— l)^ -1 )/ 2 is present because by convention Hl{t) ~ 
t 1-7 and 2-^1(2) ~ - 2;1_c i n a neighborhood of t = (resp. z — 0), where the 
principal branches are meant. The corresponding reduction of if/ to 2 -Pi is 
trivially valid (both sides are constant functions of £, and equal unity). 

Working out the number of rational substitutions z = R(t) that may be ap- 
plied to trivial Heun equations, where R is of the form M 2 o R 1 o Mi, is a 
useful exercise. There are seven subcases of Theorem 13. 11 i.e., choices for the 
polynomial R±. Each subcase allows d to be chosen from an orbit consisting 
of m cross-ratio values: m = 3 in the harmonic subcases la and lc, m — 2 
in the equianharmonic subcases 2a and 2d, and m = 6 in the others. In any 
subcase, the 4! choices for Mi are divided equally among the m values of d, 
and there are also 3! choices for M2. So each subcase yields (4!/m)3! rational 
substitutions for each value of d, but not all are distinct. 

To count distinct rational substitutions for each value of d, note the following. 
R will map t — 0, 1, d, 00 to z — 0, 1, 00. Each of the subcases of Theorem 13.11 
has a 'signature', specifying the cardinalities of the inverse images of the points 
0, 1, 00. For example, case la has signature 2; 1; 1, which means that of those 
three points, one has two preimages and the other two have one. (Order here 
is not significant.) In all, subcases la,lb,2b,2c have signature 2; 1; 1, and the 
others have signature 3; 1; 0. By inspection, the number of distinct mappings 
of t = 0, 1, d, 00 to z = 0, 1, 00 consistent with the signature 2; 1; 1 is 36, and 
the number consistent with 3; 1; is 18. 

Kuiken [18 supplies a useful list of the 36 rational substitutions arising from 
the harmonic subcase la, but states incorrectly that they are the only ratio- 
nal substitutions that may be applied to a trivial Heun equation. Actually, 
subcases la-lc and 2a-2d give rise to 36, 36, 18; 18, 36, 36, 18 rational sub- 
stitutions, respectively. By dividing by m, it follows that for each subcase, 
the number of distinct rational substitutions per value of d is 12, 6, 6; 9, 6, 6, 9. 
Of these, exactly one-third map t = to z = 0, rather than to z = 1 or z = 00, 
and consequently yield reductions of HI to 2 -Pi, or of HI to 2-Pi- So for each 
subcase, the number of such reductions per value of d is 4, 2, 2; 3, 2, 2, 3. 
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For example, the four reductions with d — — 1 that arise from the harmonic 
subcase la are 



HI (-1, 0; 0, + T )/2; t) (6.4a) 

= 2^(0, p/2; (l + 7 )/2; t 2 ) 

1?Z (-1, 0; 0, /3, 7, (1 + /3 + 7 )/2; t) (6.4b) 

= (-l)^" 1 )/ 2 ^ (0, (1 - + 7 )/2; (1 + 7 )/2; t 2 /(t 2 - 1)) 

HI (-1, 0; 0, /3, 1 - /3, 5; t) (6.4c) 

= 4-^Fl (0, (1 - 2(3 + <J)/2; 1 - # 4t/(t + l) 2 ) 

Hl(-1, 0; 0, /3, 1 -p, 5, t) (6.4d) 

= (-4)"^ 2 7\ (0, (1 - *)/2; 1 - P; -U/{t - if) 



The reduction (j6.4a|) . which is the only one of the four in which the degree-2 
rational function R is a polynomial, is simply the trivial (i.e., a = 0) case of the 
quadratic reduction ()3.8|) . The reduction (J6.4bj) was derived above as (jOj) . but 
fl6~4^|) and (l6~4d]) are new. They are related by composition with z t— > z/(z—l), 
i.e., by the involution in that interchanges z — 1 and z = oo. 

Remarkably, many rational reductions of trivial Heun equations to the hyper- 
geometric equation are not derived from the polynomial reductions of Theo- 
rem The following curious degree-4 reduction is an example. The rational 
function 

takes t — 0, 1, d = 2, oo to z — 0; and t = l±itoz = l,oo (both with 
quadruple multiplicity). By Proposition ^. 51 a trivial Heun equation with d — 2 
will be reduced by Q to a hypergeometric equation iff Q maps exponents 
to exponents. This constrains the singular points t = 0, 1, d, oo to have the 
same exponents; which by Fuchs's relation ()2.3|) is possible only if each has 
exponents 0,1/2; which must also be the exponents of z = 0. Also, since 
t — 1 ± i are ordinary points of the Heun equation, with exponents 0,1, 
the exponents of the hypergeometric equation at z — l,oo must be 0,1/4. 
It follows that on the level of solutions, the reduction is 



Hl(2, 0; 0, |, ~, |; t) = (i/4)^^ (o, 1; |; ^-i^l ' 



(6.6) 



where the normalization factor (i/4) 1 / 2 follows from the known behavior of 
the functions Hl(t) and 2-^1 (z) as t — > and 2 — > 0. 

Each of the preceding rational reductions of a trivial HI to a 2-^1 can be 
converted to a rational reduction of a nontrivial HI to a 2-F1, by using the 
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definitions (|2.9j) . (j2.11|) of Hl, 2 F\. For example, implies 

(6.7) 
8itQ-l)(t-2) \ 

(t-i + i) 4 / 

The equality ()6.7|) holds in a neighborhood of t = (both sides are real when 
t is real and sufficiently small). This reduction is not related to the previously 
derived harmonic reduction ()3.5a|) . in which d = 2 also. The pair (d,q/aj3) 
here equals (2, |), which is not listed in Theorem 13.71 

The formula (J6.7j) is a reduction of a nontrivial HI to a 2-F1, but of a more 
general type than has been considered in this paper. The underlying reduction 
of the Heun equation to the hypergeometric equation (j^J includes a linear 
change of the dependent variable, resembling a complicated F-homotopy, in 
addition to a rational change of the independent variable. 
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